Non-Markovian dynamics of a qubit 
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In this paper we investigate the non-Markovian dynamics of a qubit by comparing two generalized 
master equations with memory. In the case of a thermal bath, we derive the solution of the recently 
proposed post-Markovian master equation [A. Shabani and D.A. Lidar, Phys. Rev. A 71, 020101(R) 
(2005)] and we study the dynamics for an exponentially decaying memory kernel. We compare the 
solution of the post-Markovian master equation with the solution of the typical memory kernel 
master equation. Our results lead to a new physical interpretation of the reservoir correlation 
function and bring to light the limits of usability of master equations with memory for the system 
under consideration. 
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I. INTRODUCTION 

The dynamics of systems interacting with their sur- 
roundings is in general very complicated. Very often, 
however, the physical systems of interest are sufficiently 
isolated from their environment to allow the use of certain 
approximations such as the weak coupling approximation 
and the Markovian approximation |l(. The former one 
assumes that the interaction between the system and the 
environment is sufficiently weak, i.e. the system is quasi- 
closed. The latter one relies on the assumption that the 
characteristic times of the system are much longer than 
those of the environment, and it always assumes the va- 
lidity of the weak coupling approximation. 

Most of the results on open systems dynamics are 
based on the weak coupling and Markovian approxima- 
tions. Recent studies have shown the limits of the Marko- 
vian description of quantum computation and quantum 
error correction [J.|J.y.|o|-|(|. Moreover, nanotechnology- 
based devices using hybrid systems, e.g. combining quan- 
tum optical and solid state systems, have been investi- 
gated and seem to be very promising for future technolog- 
ical applications 0, HJ . In order to describe decoherence 
in many solid state systems non-Markovian approaches 
need often to be used 

HE miH- Finally, a compre- 
hensive and complete understanding of the interaction 
between a quantum system and its environment, not re- 
lying on the weak coupling and/or Markovian approxi- 
mations, is crucial in order to clarify fundamental issues 
such as the quantum-classical border, and in order to gain 
new insight in the dynamics of quantum systems which 
are not in thermal equilibrium. 

Outside the region of validity of the Markovian ap- 
proximation the master equation describing the dynam- 
ics cannot be usually cast in the well known Lindblad 
form [la. Il4|. This fact has several consequences: one 
for all, complete positivity of the dynamical map E3 
is not guaranteed anymore, and even positivity may be 
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violated. The latter of this properties, i.e. positivity, 
is necessary to guarantee the statistical interpretation of 
the density matrix, while the former one ensures that the 
time evolution in the system-environment total space is 
unitary . 

Non-Lindblad master equations are much more diffi- 
cult to solve, both analytically and numerically, than 
Lindblad ones, and they may lead to non-physical behav- 
iors such as, e.g., violation of positivity of the dynami- 
cal map [See, e.g., Refs. [3 El El]- 

The break down 

of the positivity condition stems from the phenomeno- 
logical nature of most of the non-Markovian approaches. 
Exact generalized master equations indeed, by definition, 
do not violate neither positivity nor complete positivity 
but they are generally far too complicated for providing 
useful means for studying the system dynamics. 

In this paper we focus on a basic model of an open 
quantum system with memory, i.e. a two-level system 
(qubit) interacting with a bosonic thermal reservoir. We 
apply a recently proposed post-Markovian approach E3 
which interpolates between the exact Kraus map and 
the Markovian dynamics. We compare the solution of 
the post-Markovian generalized master equation with the 
solution of the typical memory kernel master equation 
P, |2(| . For a specific physically interesting form of the 
reservoir spectral density it is possible to derive an exact 
solution starting from a microscopic description of the 
total system, i.e. system plus reservoir. This fact allows 
us to make a comparison between the phenomenologi- 
cal memory kernel approaches and the exact microscopic 
approach. From this comparison a new physical interpre- 
tation of the reservoir correlation function will emerge. 
Finally the limitations of the memory kernel approaches 
will be underlined and issues related to the loss of posi- 
tivity of the dynamical maps will be carefully analyzed. 

The paper is structured as follows. In Sec. II we recall 
the post-Markovian master equation and we present the 
solution for the case of a qubit interacting with a thermal 
reservoir. In Sec. Ill we recall the solution of the typi- 
cally used generalized master equation with memory, we 
compare it to the post-Markovian solution, and we an- 
alyze the "non-physical" region of the parameters space 
where positivity breaks down. In Sec. IV we consider 
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the exact solution and in Sec. V we present conclusions. 



II. POST-MARKOVIAN MASTER EQUATION 
FOR A QUBIT 

A. Interpolating the Kraus and Markovian 
dynamical maps 

Very recently a new general post-Markovian master 
equation has been presented [13. An interesting fea- 
ture of this phenomenological master equation is that, 
by construction, it interpolates between the generalized 
measurement interpretation of the exact Kraus operator 
sum map, and the continuous measurement interpreta- 
tion of the Markovian dynamics. 

It is worth reminding that the dynamics of open quan- 
tum systems may be described equivalently either by 
means of the density matrix satisfying the master equa- 
tion or by means of a stochastic wave function which 
is the solution of the stochastic Schrodinger equation 
unravelling the dynamics Q. For Markovian dynamics 
there exist a physical interpretation for the time evolu- 
tion of the stochastic wave functions (quantum trajecto- 
ries). Indeed it has been shown that the quantum tra- 
jectories describe the time evolution of the system con- 
ditioned to a continuous measurement of the environ- 
ment pll . |22| . Different types of detection schemes of the 
environment (photon counting, homodyne and hetero- 
dyne detection) correspond to different unravellings (dif- 
ferent types of stochastic Schrodinger equations) . At con- 
trast, in the case of non-Mar kovian dynamics, it has been 
shown that quantum trajectories do not have a physical 
interpretation |23j | although attempts to find an interpre- 
tation in extended Hilbert states have been performed 
[2^ | . In more detail, it turns out that the stochastic wave 
function at time t represents the state the system would 
be in at that time if a measurement was performed on 
the environment at that time, and yielded a particular 
result (generalized measurement interpretation). How- 
ever, the wavefunction at time t does not have any link 
with itself at times less then t, and therefore there cannot 
be any physical interpretation of the quantum trajecto- 
ries, for non-Markovian systems [2^| . The reason can be 
traced back to the finite correlation time characterizing 
the non-Markovian bath. 

Now, since the post-Markovian master equation actu- 
ally interpolates between the exact dynamics and the 
Markovian dynamics, an analysis of the time evolution 
according to this equation may give new insight in the dy- 
namics of non-Markovian systems, in the possible physi- 
cal phenomena taking place when the Markovian approx- 
imation fails, and in their role in the break down of the 
continuous measurement interpretation. This is in fact 
what we are going to investigate in the rest of the pa- 
per. Moreover, we will study the usefulness of the post- 
Markovian master equation presented in Ref. |19| for the 
description of open quantum systems, comparing it to 



other common non-Markovian approaches. To this aim 
we consider the basic open quantum system, e.g. a two- 
level atom or qubit, interacting linearly with a quantized 
bosonic reservoir at T temperature. 



B. Post-Markovian master equation for the qubit 

The general form of the post-Markovian master equa- 
tion introduced in Ref. [19( is the following 



dt 



dt'k(t')exp(£t')p(t - 1'), 



(1) 



where p(t) is the density matrix of the reduced system, 
fc(t') is the memory kernel, and C is the Markovian Li- 
ouvillian. For the system here considered the Markovian 
Liouvillian is given by pl| 
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with 70 the phenomenological dissipation constant, N the 
mean number of excitations of the reservoir, and a± the 
spin inversion operators. In Appendix A we recall the 
form of the solution of the Markovian master equation 
via the damping basis method. 

In the following we will focus on a widely use form of 
memory kernel, namely the exponential memory kernel 



k{t) — 76 



-7* 



(3) 



It is worth stressing that k(t), which hereafter we will 
call the Shabani-Lidar memory kernel, is a quantity 
introduced phenomenologically in Ref. |l9j |. There- 
fore, it should not be confused with the memory func- 
tion, appearing in the second order approximation of the 
Nakajima-Zwanzig equation, which is related to the spec- 
tral density of the reservoir [See, e.g., 0, p. 465]. In order 
to understand the meaning of the Shabani-Lidar mem- 
ory kernel we recall that, in the measurement scheme 
approach to open quantum systems, the post-Markovian 
master equation describes a situation in which a mea- 
surement of the environment at a time t' is followed by 
a Markovian evolution, described in terms of continuous 
measurements of the environment at times t > t' . The 
time t' characterizes the bath memory effects. In this 
picture, the Shabani-Lidar memory kernel is a function 
which assigns weights to different measurements (select- 
ing different p(t')) 0. Now, having in mind these defi- 
nitions, and remembering that that the post-Markovian 
master equation is phenomenological, one might wonder 
which is the relationship between the memory kernel of 
the post-Markovian master equation and the Nakajima- 
Zwanzig memory function (also known as correlation 
function) which is related to the reservoir spectral den- 
sity. This question will be addressed in Sec. IV. 
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C. Analytic solution 

By applying the method described in [l^ , and recalled 
in Appendix B, to a two- level system whose ground and 
excited states are |1) and |2), respectively, we derive the 
following solution of the post-Markovian master equation 



Pit) 



£(R,t) U1-P22 



1 



2JV + 1 



1 



2N + 1 



[p\2<y+ +P21V-] > ,(4) 



where /Oy — (i\p(0)\j), with i,j = 1,2. In the previous 
equation 



and 
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r(R) 
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(6) 



with R = IA2I/7, the eigenvalue A2 being the one derived 
for the Markovian master equation [see Eq. l|28fl]. and 
t = yt. 

Equation© is valid only for r(R) < 1 and r(2R) < 1. 
When r(R) > 1 and/or r(2R) > 1, the form of the time 
dependent coefficients £(i2,r) and/or £(2R, r) appearing 
in Eq. Q is obtained from Eq. (JSJby substituting sinh[.] 
and cosh[.] with sin[.] and cos[.], respectively. 

Let us focus on the case of a zero-temperature reser- 
voir. In this case the solution of the post-Markovian mas- 
ter equation takes the form 



p{r) = \{I+[2P 1 {r)-l]a z 



+ [2p 12 (T)a++h.a.}}, (7) 



where the time evolution of the ground state population 
Pi(r) = pix(r) and of the coherences Pi2(t) = p^x (r) is 
given, respectively, by: 



Pi(r) = PiiZ(R,t), 
PM = ^p l2 £(2R,T) 



(8) 
(9) 



We note that, for the zero-temperature case here consid- 
ered, A2 = —70, and therefore R = 70/7- 



D. Markovian limit 

We conclude this section by showing the Markovian 
limit of the post-Markovian solution. To this aim we 
firstly notice that Eq. JSJ can be recast in the following 
simplified form 



£(R,t) 



1-R ' 



(10) 



For times f > I/7 = tr (coarse graining in time) and 
for R <C 1, i.e. r = I/70 3> tr, the previous equations 
become £(R,t) ~ e"^ 2 ' 4 , while £(2i?,i) ~ e~ 2|A2|t , and 
one reobtains the Markovian dynamics. The approxima- 
tion ro = I/70 ^ tr amounts at saying that tr is much 
smaller than the characteristic time of the system ro . In 
the following we will call tr the characteristic time of the 
reservoir. 



III. GENERALIZED MASTER EQUATION 
WITH MEMORY 



A. Phenomenological master equation 

Let us now consider the typical phenomenological mas- 
ter equation with memory kernel having the form 



(5) 



% = I dt'k{t')C P {t - t'). 

"it Jo 



(11) 



This generalized master equation takes into account the 
previous "history" (0 < t' < t) of the density matrix p(t) 
by means of the phenomenological memory kernel k(t'). 
Differently from time-convolutionless approaches [l| , this 
leads to a master equation which is not local in time. For 
specific forms of the memory kernel, as the exponential 
one considered in this paper, this master equation can be 
solved by means of the Laplace transforms. 

The memory kernel master equation given by Eq. 
(fllf) has been used in Ref. 01 f° r studying the non- 
Markovian dynamics of a quantum harmonic oscillator 
interacting with the vacuum. There the authors have 
found that, for an exponential memory kernel such as 
the one given by Eq. ©, the positivity of the density 
matrix is violated for certain values of the phenomenolog- 
ical decay constants. In Ref. |^ the generalized master 
equation with memory given in Eq. (|llf> has been ana- 
lyzed for a two-level atom in the presence of telegraphic 
noise, and conditions for complete positivity have been 
presented. In the present paper we consider this model 
for the two-level atom interacting with a bosonic reser- 
voir at T temperature, focussing in particular, for the 
sake of simplicity, on the zero temperature reservoir. 

Similarly to the case of the post-Markovian master 
equation, one can solve Eq. (|1 1|) by taking its Laplace 
transform, determining the poles, and inverting the solu- 
tion in the standard way. Using the damping basis given 
by Eq. f2T| . with the eigenvalues given by Eq. 
we find that the solution has the same form of Eq. J3J, 
with the only difference that the quantities t) is now 
given by 



t(R,r) 



■. sinh 



+ cosh 



;v1T 



•4R| 
^lR\\} 



V|l-4i?| 



(12) 
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By comparing Eq. (JSJ and Eq. I|12f> . with the help of 
Eq. ©, one easily sees that this amounts at assuming 
R — ja/j^ 1 in Eq. © . Stated another way, as one can 
see directly form Eqs. and the memory kernel 
master equation is a special case of the post-Markovian 
master equation, in the limit in which the system charac- 
teristic time Tn is much bigger than the reservoir correla- 
tion time tr [l!j . We note that, for 4R > 1 and 8R > 1 
the form of the time dependent coefficients £(R, r) and/or 
£ (2i?, r) , respectively, is obtained from Eq. (fT2f by sub- 
stituting sinh[.] and cosh[.] with sin[.] and cos[.]. 

We conclude this section noting that, as for the Marko- 
vian dynamics, the solution of the master equation with 
memory, given by Eq. can be obtained from the 

solution of the post Markovian master equation in the 
limit R = 70/7 *C 1. However, contrarily to the Marko- 
vian dynamics, no coarse graining in time has been made, 
and therefore the solution of Eq. describes correctly 
the short time dynamics, t <C tr, characterized by non- 
negligible system-reservoir correlations. 

B. Positivity of the dynamical map 

As we have already mentioned in the Introduction, 
since the master equations given by Eqs. ||TJ and 
are not of Lindblad-type, the positivity of their corre- 
sponding dynamical maps is not guaranteed. The break 
down of the positivity condition means that the density 
matrix loses its statistical interpretation, its eigenvalues 
becoming negative. This is hence a sign of failure of the 
phcnomcnological master equations with memory. 

An analysis of the positivity condition for a master 
equation of the form of Ea. i|llJ) has been presented, in 
the case of a damped harmonic oscillator, in Ref. jl7| . 
There the authors have found that positivity is always 
violated for sufficiently high values of the phenomeno- 
logical decay constant. The positivity condition of the 
memory kernel master equation has been also studied in 
Ref. [26j for different types of memory kernels, includ- 
ing the exponential one, generalizing the results of |l7j . 
It is therefore not surprising that, as we will see in the 
following, we obtain the same result of Ref. [I3 for the 
same form of master equation, with an exponential mem- 
ory kernel, when the system is a qubit. As we will show, 
however, the post-Markovian master equation exhibits a 
strikingly different behavior in respect of the positivity 
condition. 

In order to study in more detail the positivity condition 
for both the two memory kernel master equations consid- 
ered in this paper, it is more convenient to rewrite the 
solutions in terms of the Bloch vector w — {w x , w v , w z }. 
The qubit density operator at time r = jt, given by 
Eqs.QJ and Eq. (JjJJ in the post-Markovian case, and by 
Eqs.QJ and (|12[1 in the memory kernel case, can be recast 
in the form 

p(r) = ^[I + w(r)-ff], (13) 



with <t = {ax-, c y , Cz}, and 

w x (t) = £(2R,T)Re{p 12 ]=Z(2R,T)w x (0) (14) 
w v {t) = -Z(2R,T)Im[p u ]=Z(2R,T)Wy(0) (15) 
w z (t) = 2P 1 Z(R,t)-1 = P 1 (t)-P 2 {t), (16) 

where P\{t) is given by Eq.(JHJ), P2O") = 1 — Pi( T ): an d 
£(i?, r) is given by Eq. 112JI (post-Markovian) or Eq. 
lj5Jl (memory kernel) . The dynamical map describing the 
evolution of the qubit is positive if and only if the density 
operator evolves only to states inside or on the Bloch 
sphere. Therefore the positivity condition, in terms of 
the Bloch vector components, simply reads |mj(r)| < 1, 
i = x.y, z. By looking at Eqs. (JTIJ-JlU) and Eq.©, 
one sees immediately that the conditions Wi(f) < 1 are 
equivalent to |£(-R,t)| < 1 and \£(2R,t)\ < 1. 

Let us begin considering the post-Markovian master 
equation. It is straightforward to prove that £(i?, r) is 
a positive and monotonically decreasing function of R 1 
therefore, it is sufficient to investigate the conditions for 
which \£(R,t)\ < 1, since \£(2R,t)\ < \£(R,t)\ for all 
values of R and r. As shown in Appendix C, it turns 
out that the post-Markovian dynamical map never vi- 
olates positivity. On the contrary, for the case of the 
memory kernel master equation given by Eq. (|llf> . a 
close look to Eq. ljT2*|) tells us that \£(R, t)\ < 1 only for 
AR < 1, i.e. 470 < 7. We remind that one may derive 
the memory kernel master equation given by Eq. Hll|) 
from the post-Markovian master equation given by Eq. 
(|T|) in the limit 70 < 7. Having this is mind it is not 
surprising that the positivity condition breaks down for 
sufficiently high values of 70, for the system here con- 
sidered. The study of positivity, therefore, suggests that 
the post-Markovian master equation is somehow "more 
fundamental" than the memory kernel master equation. 

IV. EXACT DYNAMICS 

In this section we will analyze two different aspects 
characterizing the non-Markovian dynamics of a quibit 
described by means of memory kernel master equa- 
tions. The first aspect stems from the microscopic deriva- 
tion, using the Nakajima-Zwanzig formalism, of the phe- 
nomenological memory kernel master equation given by 
Eq. I|ll|) . The microscopic derivation will allow us to 
link the Shabani-Lidar memory kernel to the correlation 
function and to the reservoir spectral density, in the limit 
7o < 7- 

The second aspect described in this section is related 
to a physically interesting specific microscopic model of 
non-Markovian dynamics of a qubit for which an exact 
solution does exist. The existence of an exact solution 
allows to make a comparison with the predictions of the 
post-Markovian and memory kernel solutions, and hence 
to study the limits of both these approaches. Moreover, 
keeping in mind that the post-Markovian approach can 
be seen as an interpolation between the exact Kraus map 
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and the Markovian one, one may gain new insight in the 
reason why non-Markovian quantum trajectories lack of 
a physical meaning. 



A. Microscopic derivation 

Let us begin with the microscopic derivation of the 
memory kernel master equation given by Eq. Ijlll) . Wc 
consider the following microscopic Hamiltonian for the 
total system, i.e. system plus reservoir, 



with 



H = H + H T , 



H = + y^fcbfcbfc 

fc 

Hi = a+B + a-B^, 



(17) 



(18) 



where B = X^fcfffc^fc- The transition frequency of the 
two- level system is denoted by w , the index k labels 
the different modes of the reservoir with frequencies u>k, 
b\, and indicate the creation and annihilation oper- 
ators, and are the coupling constants. Following 
the typical approach for the derivation of master equa- 
tions for the reduced density matrix we firstly write the 
von Neumann master equation for the total system in 
the interaction picture, then we trace over the environ- 
mental degrees of freedom, under the assumptions that 
Tr,R \H[{t), pt(0)] = 0, with pt the density matrix of 
the total system, and pr(0) = p(0) (g> pr(0), with Pr(0) 
the density matrix of the reservoir. To the second order 
in perturbation theory (Born approximation), we obtain 
the following integro-differential equation for the reduced 
density matrix 



dp(t) 
dt 



dt'TiR [Hj(t) 7 [Hj(t'), p(t') ® p R }} . (19) 



Having in mind Ea. l|18[) it is then straightforward to re- 
cast Eq. i|19|) in the same form of Eq. <|ll|l with 



k(t') = Tr R [B(0)B\t') P R] 



(20) 



In the previous equation B(t) — J2k 9^>k exp(— iu>kt) is 
the reservoir operator appearing in Eq. I|18|) . in the in- 
teraction picture, and J(oj) = J2k9k3( UJ -Wfc)/(2mi : w/ ! ), 
with rrik masses of the oscillators of the reservoir. From 
the previous definition of k{t') one sees clearly that this 
function describes temporal correlations of the reservoir 
operators B, and therefore it is commonly known as the 
reservoir correlation function. In the second line of Eq. 
(|20|l we have introduced the so-called spectral density of 
the reservoir J(uj), which is therefore simply the Fourier 
transform of the correlation function. An exponential 
correlation function of the form of Eq. J2J corresponds to 



a Lorentzian spectral density as the one typical of cavity 
quantum electrodynamics, for a reservoir in the vacuum 
state. 

In the rest of this section we will focus on the following 
specific physical system: a two level atom interacting 
resonantly with a quantized mode of an empty high Q 
cavity. Assuming that the two-level atom is in resonance 
with the cavity mode, the reservoir spectral density is the 
following 



1 



J(u) = - 



ToA 2 



7T (lUq — Uj) 2 + X 2 



Using Eq. (20J), we get 



k(t') = 7 Ae 



and the master equation (|19fl becomes 



dp 
~dl 



= / dt'Ck(t')p(t-t'), 



(21) 



(22) 



(23) 



with k{t') given by Eq. (|2*2)l . and 7o£ = (70/70)^ , where 
C is given by Eq. @, with N = 0. A direct compari- 
son with Eqs. (J3J) and i|ll|) clearly shows that the master 
equation derived using second order perturbation theory 
starting from the microscopic description above coincides 
with the phenomenological memory kernel master equa- 
tion given by Ea. l|ll|l . provided that A = 7, 70 = 70. 



B. Physical interpretation of the system and 
reservoir parameters 

The microscopic derivation allows to give a physical 
interpretation of the decay constants 70 and 7. Indeed, 
Eq. (|21|l tells us that 70 measures the strength of cou- 
pling between the two-level atom and the vacuum reser- 
voir, and hence the system characteristic time tq is de- 
termined only by the system-reservoir coupling strength. 
The reservoir correlation time tr = I/7 is simply given 
by the inverse of the width of the Lorentzian spectral 
density. 

Having this in mind, the reason of the violation of pos- 
itivity for the master equation given by Eq. i|ll|) becomes 
clear. Indeed, such master equation correctly describes 
the dynamics of the system only when second order per- 
turbation theory is valid, i.e. for weak system-reservoir 
coupling. The conditions in correspondence of which pos- 
itivity is violated, i.e. when 4i? > 1 and 8R > 1, with 
R = 70/7, however, correspond to strong couplings be- 
tween system and reservoir. 

IhTaI 



Finally, we remind that, as noticed in Sec. 
Eci. (|19|) is the limit for 70/7 <C 1 of the post-Markovian 
master equation given by Eq.([Q|. As a consequence, in 
this limit, the Shabani-Lidar phenomenological memory 
kernel k(t') coincides with the Fourier transform of the 
reservoir spectral density, as given by Eq. (|20|l . This fact 
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allows to give a new physical interpretation to the cor- 
relation function of the reservoir in terms of generalized 
measurement by the environment. Indeed, it turns out 
that the reservoir correlation function acts as a weighting 
time distribution function, assigning weights to different 
measurements selecting different p(t'). 



C. Exact solution 

The physical system we consider in this section is one 
of the few open quantum systems amenable of an exact 
solution [22]] . For a spectral density of the type given by 
Eq. (|21|) , the exact density matrix has the following form 



1] 



Pit) 



Pi(r) Pi2(r) 
P* 12 (t) l-Px(r) 



(24) 



with 



Pi(r) 



Pi(0)e- T |cosh 
1 



V|l-2i?| 



: sinh 



V\l^2R\l\}\ (25) 



for 2R < 1. If 2R > 1, Pi(t) has the same form of 
Eq. 125fl provided that the substitutions cosh[.] — > cos[.] 
and sinh[.] — > sin[.] are made. In Fig^we compare the 
time evolution of the exited state population as predicted 
by both the post-Mar kovian master equation [see Eqs. 
© and (jSJ)] and by the memory kernel master equation 
[see Eqs. (JSJl and l|12fl ] with the exact dynamics [see Eq. 
iPSJ], in correspondence to three different values of the 
parameter R — 70/7. We assume that the initial state of 
the two-level system is the exited state. From the figure 
one can see that while the memory kernel approximated 
dynamics does violate positivity for strong enough cou- 
plings (R = 5 and R = 1), the post-Markovian dynamics 
is always positive. However, both the post-Markovian 
solution and the second-order solution approximate well 
at all times r the exact dynamics only for small values 
of R, e.g. for small couplings. Therefore, the specific 
example here considered shows that there exist situa- 
tions for which there is actually no advantage in using 
the post-Markovian approach when compared to second 
order perturbation theory. The reason why the post- 
Markovian approach fails in describing the dynamics of 
the system even for intermediate couplings is related to 
the way in which such a master equation is derived. Let 
us recall the physical meaning of the post-Markovian ap- 
proach in terms of generalized measurement interpreta- 
tion. The derivation of the post-Markovian master equa- 
tion assumes that, after the time t' at which the gen- 
eralized measurement by the environment is performed, 
the evolution of the system is Markovian. The time dis- 
tribution of the instants t' at which the measurement 
is performed is given by the memory kernel. Therefore 
one should expect such master equation to be valid for 
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FIG. 1: Time evolution of the excited state probability for 
three different values of R, i.e. R = 5, R = 1, R = 0.05. 
The dimensionless time is r = 'yt. The solid line indicates 
the dynamics for the post-Markovian master equation, the 
dashed lines indicates the exact dynamics and the dotted line 
the dynamics of the memory kernel master equation derived 
using second order perturbation theory. 



To > Tr. Indeed, for To < tr the assumption that the 
dynamics after t' is Markovian would not be well justified 
since the reservoir correlation time would then be longer 
than the Markovian dissipation time, and this would in- 
evitably lead to a non-negligible feedback of the environ- 
ment to the system. We see from Eq. (|25f) that already 
for R — 70/7 > 0.5 the exact dynamics of Pi(t) shows 
an oscillatory behavior. These oscillations may be seen, 
in a completely quantum approach, as virtual absorption 
and re-emission of the same quantum of energy from the 
environment. The description of these quantum phenom- 
ena cannot be present in the post-Markovian approach. 
It is exactly the appearance of these virtual exchanges of 
energy which does not allow to give a physical interpreta- 
tion to the single trajectories for strongly non-Markovian 
systems, since there seems to be no way for a single phys- 
ical trajectory to describe a virtual process. 



V. SUMMARY AND CONCLUSIONS 

In this paper we have taken into consideration two 
models of generalized master equations with memory, 
and we have applied them to the description of the non- 
Markovian dynamics of a qubit interacting with a quan- 
tized bosonic reservoir in thermal equilibrium. For the 
case of an exponential memory kernel we have compared 
the solution of the recently proposed post-Markovian 
master equation with the solution of the typical master 
equation with memory kernel. We have demonstrated 
that, for the system considered, the post-Markovian ap- 
proach never violates positivity, contrarily to the mem- 
ory kernel master equation. We have then seen that the 
memory kernel master equation coincides with the second 
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order expansion of the exact Nakajima-Zwanzig general- 
ized master equation. Since the memory kernel master 
equation is the limit of the post-Markovian master equa- 
tion for 70 -C 7, it is possible to give a generalized mea- 
surement interpretation to the correlation function of the 
reservoir. Finally we have considered the following phys- 
ical implementation of the system: the qubit describes 
the excited and ground electronic state of an effective 
two-level atom crossing a high Q cavity; the reservoir 
is formed by the quantized modes of the high Q cavity 
which are distributed according to a Lorentzian peaked 
at the atomic Bohr frequency. This physical system, typ- 
ical of cavity QED, represents one of the few examples 
of exactly solvable open quantum systems. The compari- 
son between the exact dynamics and the post-Markovian 
and memory kernel solutions shows that there exist situ- 
ations in which the post-Markovian approach does not 
present any advantage over the second order approxi- 
mated memory kernel master equation. The reason is 
traceable back to the fact that, by derivation, the post- 
Markovian master equation cannot describe accurately 
situations for which the characteristic time of the reser- 
voir tr is greater than the characteristic time of the sys- 
tem t . When tr > t the dynamics is characterized 
by virtual exchanges of energies between the system and 
the environment which cannot be described by the post- 
Markovian approach. Such virtual processes, absent in 
the Markovian dynamics, seem to be responsible for the 
lack of a physical interpretation of single quantum tra- 
jectories in terms of continuous measurements performed 
by the environment. 
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Appendix A 

The damping basis method allows to solve the mas- 
ter equation given by Eq.J5J by solving the eigenvalue 
equation 

Cpx = X Px . (26) 

It turns out that the damping basis is |2^| 

Pa, : {p\ x = vo,P\ 2 = °"z,Pa 3 = o + ,Pa 4 = f-}, (27) 

with (To = \[I — o' z /(2N + 1)], a± — u x ± i<jy, and 
o x ,<jy,cr z the Pauli matrices. The corresponding eigen- 
values are 

\ : {Ax = 0, A 2 = -2 7o (iV+l/2), A 3 = A 4 = - 7o (jV+l/2)}. 

(28) 



The density matrix p(t) can then be written, in the 
damping basis, as follows 



(29) 



with c\ i = Tr {paiP(O)}, where p\ i is the dual damping 
basis. Inserting the values of Eqs. H27|) - (|28() into Eq. 
(|2"§|) . one gets 



Pit) = - 2 {I- 
1 



1 



2N+1 



2N + 1 



e -7o(JV+l/2)t 
Oz H -Z [P12C+ + P2lCT-|:}0.) 



Appendix B 

In this Appendix we recall the main steps to derive the 
general solution of Eq. Q , as demonstrated in , and 
we carry out the derivation for the case of a qubit in- 
teracting with a quantized thermal reservoir. The initial 
step to solve the post-Markovian master equation is the 
derivation of the damping basis for the Markovian case 
(see Appendix A). As in the previous Appendix, we de- 
note with {Xi} the complex eigenvalues and with {/?>;} 
and {p\i } the damping basis and its dual, respectively. 
Then we write the density matrix as follows 

p(t) =p Ji (t)p Xi . (31) 

Taking the Laplace transform of Eq. JTJ one gets 

C 



sp(s) - p(0) 



k(s) 



-C 



P(s) 



(32) 



where * denotes the convolution. Taking the Laplace 
transform of Eq. (|31|l and using the previous equation 
one obtains 



sp,i(s) - Hi(0) = Xik(s - Xi)jli(s) 
and transforming back 
1 



(33) 



Pi = Lap 1 



Mi (0)=&(i) M i(0). (34) 



.s - Xik(s - Xi)_ 

The coefficients £j may be calculated once fixed C and 
k(t) (see ^3 )> therefore one gets 



p(t) = ^i(t) f x i (0) P X i 



(35) 



For the case of a qubit interacting with a T- 
temperature reservoir, the damping basis is given by Eq. 
(|27|l . Assuming an exponential memory kernel, and using 
Eq. (|2"5|> . we have solved Eq. (|3~4"}l obtaining 



aw = 1, 



(36) 



with ^(R^t) given by Eq. (|3J). Inserting the previous 
equations into Eq. l|35f) . one gets Eq. Q. 
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Appendix C 

In this Appendix we demonstrate that the post- 
Markovian master equation for a qubit never violates the 
positivity condition for an exponential memory kernel. 
We have seen in Sec. IIII Bl that the positivity condition 
amounts at < £(R, r) < 1. 

Let us firstly show that £(IZ, r) < 0. By looking at Eq. 
l(TU)l. and remembering that for the zero temperature case 
I ^2 1 = 7o j one sees immediately that this corresponds to 
prove that 



The first set of inequalities is always satisfied since when 
R < 1 then e (1 ~ fl ) T > R at all times r. Similarly the 
second set of inequalities is always satisfied since when 
R > 1 then e - ( H-1 ) T < 1 < R at all times r. 



We now show that r) 
£(R, r) > we have 



-Rt 



Re 



l-R 



< 



< 1. Since R > and 



Re 



1 - R 



= e~ T < 1.(38) 



R > 0, 
i? < 0, 



e- RT - Rer T > 0; 
e~ RT - Rer T < 0. 



(37) 
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